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Abstract 

We solve the general one-dimensional Dirac equation using a "Poincare Map" approach which 
avoids any approximation to the spacial derivatives and reduces the problem to a simple recursive 
relation which is very practical from the numerical implementation point of view. To test the 
efficiency and rapid convergence of this approach we apply it to a vector coupling Woods-Saxon 
potential, which is exactly solvable. Comparison with available analytical results is impressive and 
hence validates the accuracy and efficiency of this method. 
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The Dirac equation describes a relativistic particle of spin one-half at high velocities (below the 
threshold of pair production) [TJ. It describes the state of electrons in a way consistent with special 
relativity, requiring that electrons have spin 1/2, and predicting the existence of an antiparticle partner 
to the electron (the positron). The physics and mathematics of the Dirac equation is very rich, it is a 
first order matrix linear differential equation whose solution is a 4-component wavefunction (a spinor) . 
Nevertheless, it was hard to find nontrivial exact solutions of this equation. Until 1989, there was 
only one nontrivial exact solution of the Dirac equation for the very important Coulomb problem. By 
convention, a nontrivial exact solution of the Dirac equation is a solution that, in the non-relativistic 
limit, reproduces a known solution of Schrodinger equation with that specific potential and hence 
carries the same relativistic potential name. 

In this regard, Moshinsky and Szczepaniak [2] in 1989 were able to formulate and solve the relativis- 
tic oscillator problem (Dirac-Oscillator). During the last decade, Alhaidari has introduced an effective 
approach to the solution of the Dirac equation for spherically symmetric potentials [31 El El El E]- 
His method was initiated by the observation that different potentials can be grouped into classes; 
for example, the non-relativistic Coulomb, oscillator and S-wave Morse problems constitute one such 
class. Therefore, the solution of two problems in one class implies solution for the remaining one. By 
this method, the S-wave Dirac- Morse problem was formulated and solved [3J. 

Using the above-mentioned method, other potentials were treated; among these are Dirac-Scarf, 
Dirac- Rosen- Morse I & II, Dirac- Poschl- Teller, Dirac-Eckart [3], Dirac-Hulthen, and Dirac- Woods- 
Saxon potentials [8]. On the other hand alternative methods such as quasi-exactly solvable problems at 
rest mass energies with power-law relativistic potentials were investigated following the same procedure 
[5]. Orthogonal polynomials were also used to find series solutions of Dirac equation for scattering and 
bound states [6]. In this work we will be interested in solving the one-dimensional (ID) Dirac equation 
in general using a Poincare map which enable us to treat exactly the spacial derivative operator while 
the only approximation is made upon discretizing the scattering potential. The approach is very 
efficient and converges rapidly. To exhibit the efficiency of our method we use it to solve the vector 
coupling Woods-Saxon potential and compare our results for the transmission coefficient with the 
analytic one. The accuracy and easy implementation of our method is impressive which suggest its 
suitability in dealing with ID Dirac equation for an arbitrary potential. 

We consider a one dimensional Dirac equation for a particle of mass m, subject to a vector potential 
coupling V(x). In unit system H = c = 1, this particle is governed by the Dirac Hamiltonian 



with a x and a z being the Pauli matrices. The stationary eigenvalue equation can be written explicitly 
in the following form 



where tp(x) = (ip 1 ip 2 ) 1 is the eigenspinors of two components. Generally speaking for an arbitrary 
potential V(x), this equation is difficult to solve analytically except for a very limited number of 
solvable potential as classified by super symmetric quantum mechanics jllj . Alternatively, in this work 
an iterative method that was used a lot for ID Schrodinger equation, the Poincare Map approach |12j . 
will be adapted to our relativistic problem. Using this approach we are going to show that the above 
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ID Dirac equation can be replaced by a Poinccare map associated with the above wave equation. This 
approach will enable us to solve the Dirac equation and generate the spinor wavefunction iteratively, 
a method very suitable for numerical implementations. In addition, we will be able to use the transfer 
matrix approach to compute the transmission coefficient. 

We consider our one-dimensional system where the particle is moving under the action of a scat- 
tering potential V(x) bound to a region of size L of our system. So our space can be decomposed 
into three major regions: the extremes are free like regions where the potential is almost zero and an 
intermediate region where the potential is V(x). Now we proceed in subdividing the potential interval 
L into N + 1 equal regions, in every region we approximate the potential by a constant local value 
V n = V(x n ) where x n = nh and h = j/^j, with N is the maximum value of n. Hence, the Hamiltonian 
in each region (n), defined by h(n — 1) < x < hn, is given by 



H n = p x a x + ma z + V n . 



(3) 



Actually we can even generalize this approach to treat spatially dependent masses in which case 
m n = m(x n ) in the above equation, however from now on we limit ourselves to a uniform mass 
distribution. In each spacial region we have a ID Dirac fermion subject to a constant potential, hence 
the solution is easily generated. The complete solution for © satisfying eigenequation with spinor 
Tp n (x) = (ipn V'n)* i n n-th space region where the potential is approximated by it local value V n 
will depend on the energy range and can be written as 



if} n {x) = A 

ip n (x) = A 
where we have set 




-ip n x 



-ip n X 



£ <V n — m 
V n - 171 < £ < V n 

V n < e < V n + m 
£>V n + m 



(4) 
(5) 
(6) 
(7) 



\s -Vn.\ — m 



\e — V n \+m 



y/\(m + V n -e)(m-V n + e)\. 



(8) 



The coefficients A n and B n are two normalization constants. The above solutions can be written in a 
compact matrix form, such as 



^„(X) = M*(X) 



A,, 



(9) 



where the matrix M*(x) is given by 



/ 



M*{x) 



i 2 a r 



PnX 



% 2 a r 



-i 2 a r 



-i 2 a. 



(10) 
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= sign(e — V n ) and s n = sign(je — V n \ —m) are the usual sign functions, equal to ± for a positive and 
negative argument, respectively. In order to treat scattering problems, and for ease of implementation, 
we consider an incident wave propagating from right to left, this amounts to change i into —i in our 
previous spinor solution. 




L n N+1 

Figure 1: Space discretization of the scattering potential V(x) 
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Figure 2: Solution space for two consecutive regions 



To obtain the desired Poincare Map we use the continuity of the spinor wavefunction at the junction 
x = x n separating the n-th and (n + l)-th region. This gives the relationship 



M n {x n ) 



L n +i(x r , 



\ B n J y B n+ i J 

We also have a relationship between M n+ \ and L n +\, which reads 



L n+1 {x n ) = M n+1 (x n ). 



(11) 



M n+ i(x n+ i) — L n+ x{x n )S n +i, S n 



s "+i) /2 Pn+1 h 



+i 



/ 1+s ^) /2 Pn+1 h 
Using the above results we can write the desired Poincare Map as 

1p n+1 (x n+1 ) = T n 1p n (x n ), T n = M n+1 (x n+ l)S n+1 M~l 1 (x n+ i) = 



(12) 



Til Tl2 
T21 T 2 2 



where the explicit matrix elements of T n are given by 



Til 



a(bc — d) 
c — d 



T12 



q(l - b) 
c — d 



T21 



a(b - l)c(a n+ i) s ^ 



and 



a = e 



^ 1+ ^) /2 a n+ i) S;+1 , 



c — d 



T22 



a(c — &d) 
c — d 



;(i+< +1 )/2 



«n+l 



b n+l 



(13) 

(14) 

(15) 
(16) 
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Now we proceed to determine the transmission amplitude t using the above recursive relationship, 
which connects A n and ip n . For x < where Vq = (i.e. left region), we can use for our transmitted 
spinor evaluated at n = 0, which when suitably normalized can be written as 




(17) 



On the other side, for x > L were Vo = (i.e. right region), we have both incident and reflected spinor 
waves. Just outside the potential region on the right hand side in the (N + 2)-th region the spinor 
wave is given by 

( i A _ ( i 



if> R (x) = A 



N+2 



-ip x 



+ B 



N+2 



\ — iao l \ icto l 

From the above notation we can easily define the transmission amplitude as 



t 



1 



(18) 



(19) 



An+2 

Hence to evaluate the transmission amplitude all we need is to find ^4jv+2 using the above recur- 
sive scheme. Our strategy now is to express Ajy + 2 in terms of the two end points spinors tpN+i = 
i/>N+l(%N+l) and ipN+2 = ?Pn+2(xn+2)- This can be easily done using our previous relations and leads 
to the form 

p i( 1+s o) /2 po h(N+2) 

(20) 



A 



N+2 



H + 1) 

i(l+«§)2 



O. 



e 2i( 1+s o) /2 p h 
(l- S l/2) 



;(i+»?)/2 



n 



(l+«S)/2 



lpN+2 



4)/- 



p h 



Summing up, we iterate the Poincae Map (|13j) to obtain the end point spinors, ipN+i and ipN+2, in 
terms of the transmitted signor ipQ. These spinors will then be injected in equations (|20p and (|19|) to 
determine the transmission amplitude. 

To test the validity of our previous approach we will now implement it for a test potential, the 
Woods-Saxon potential which has an exact analytical solution \l'6\ I14j. This potential is defined by 



V(x) = V 



6{-x) 



+ 



9{x) 



(21) 



\ _)_ g— a(x+L) ^ _|_ e a(x-L) 

where Vo is a real and positive for a barrier or negative for a well, a and L are real and positive. 
6(x) is the Heaviside step function. The analytical solutions provided in reference [T5] is used to 
evaluate the exact transmission coefficient for the potential (|2ip from the asymptotic behavior of the 
wavefunctions. This analysis leads to 



1 



E + k 
E-k 



(22) 



where A and B are given by 

r(i - 2//)r(-2z/) 



.4 



B 



T(-n -v- A)r(l - n - v + A) 
r(l - 2/j,)T(2v) 

r(- / u + v- A)r(i - /i + + A)' 



-27T/J 



-27T/J 



+ A 



+ A 



r(l + 2n)V(-2v) 



Y(n - v - A)r(l + /i - v + A) 

r(i + 2/i)r(2t/) 
r()u + i/- A)r(i + ^ + z/ + a) 1 



»c 



(23) 
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with the ratio 

D2 = r(2 j u)r(i -ji-v- A)r(-/i -v + \) 2i ^ iaLfl 

D 1 r(-2//)r(l + fi - v - X)T(fi - v + A) 

and the abbreviations fx = ^-,u = = p 2 = (E — Vo) 2 — m2 , k 2 = E 2 — m 2 have been used. 

For scattering states, of interest to us, \E\ > m ensures that k is real while the momentum p is 
real for m < E < Vo — m (the Klein range) and E > Vq + m, it is imaginary for Vq — m < E < Vq + m. 
The potential strength Vq is real and positive in our computations. In figure 3 below we show the 
shape of Woods-Saxon potential for the parameters L = 10, a = 5 and Vq = 1.2, the vertical lines 
represent the discretization of this potential. We shifted the potential to the right by an amount L 
for convenience, such a translation does not affect the physics of the problem. 




Figure 3: The Woods-Saxon potential curve V(x) shifted to the right by an amount L. 

Using the above mentioned numerical procedure (Poincare Map) we evaluate the transmission coef- 
ficient associated with the potential characterized by the above mentioned parameters (L, a, Vb) and 
a mass m = 0.4 in our atomic units. Moreover, using the same parameters as in the literature |13j . 
we have plotted the transmission coefficient in the Klein zone for a constant mass m. As shown in 
Figure 4 the solid lines correspond to the exact transmission |13j . and the dashed lines are generated 
by our Poincare iterative map. We notice from this figure, as expected, that as N increases, the dotted 
line curve converge to the exact solid line curve. We see the satisfactory agreement between our fast 
converging numerical approach and the analytical results [131 E] • 




(a) (b) 
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Figure 4: The transmission coefficient as a function of energy for two different iterations, (a): N = 600 
and (b): N = 1000 



In Figure 5 we show the transmission coefficient as a function of the potential strength Vo for two val- 
ues of the iteration parameter N. The agreement between our numerical approach and the analytical 
results is impressive. Thus we expect that we can apply our iterative approach to very general poten- 
tial, which do not lend themselves to analytical solutions. In summary, we believe that the Poincare 
Map exposed in this work is very simple and enable us to solve the ID Dirac equation for any arbi- 
trary short range potential with high accuracy and simple computational means. This approach can 
be easily extended to handle two-dimensional Dirac equation, which is playing an important role in 
describing the recently discovered graphene system [15] . 



Figure 5: The transmission coefficient as a function of Vq for two different iterations, (a): N = 400 and 
(b): N = 1000. 
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